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Swimming and pumping at low Reynolds num- tical trapping potential introduces a new control 
cd | bers are subject to the "Scallop theorem" , which parameter. We show that this system is able to 
states that there is no net fluid flow for time generate flow in a controlled fashion, character- 
reversible motions. Microscale organisms such izing the model experimentally and numerically. 
as bacteria and cells are subject to this con- 
straint, and so are existing and future artificial 

"nano-bots" or microfluidic pumps. We study i Introduction 
a very simple mechanism to induce fluid pump- 
ing, based on the forced motion of three colloidal T . , , . , . t ■ , 
. . . Insight into swimming of microorganisms and 
beads through a cycle that breaks time-reversal u , . u . j , , , . ,, 

X \ bacteria can be gained by studying the mo- 
symmetry. Optical tweezers are used to vary the ,. , , ,, , t 
J J ^ J tion at low Reynolds numbers of expenmen- 

inter-bead distance. This model is inspired by , i j , , , • , -, i , m roi n 

v J tal and theoretical model systems \UU\- uom- 

a theoretical swimmer proposed by Naiafi and -, , , u • , -, , -i • { 

^ ^ j j pared to the macroscopic world, a striking tea- 

Golestanian Phys.Rev. E, 69, 062901, 2004, , t , . ,,, F . , ' , 

1 J ' ' ' J ' ture of propulsion of these micron-scale objects, 

but in this work the relative softness of the op- t -,. , ■ ,. . , . „ ,- - u , . , , , 

r tor which inertia is typically negligible, is that 



a body striving to move has to change its shape 
with time in a non- reciprocal fashion [3] . For ex- 
ample, swimming of a magnetically driven semi- 
flexible artificial filament was recently demon- 
strated, by a wave-like motion similar to flag- 
ella [1]. In this work we realize a much simpler 
system which uses only two effective degrees of 
freedom, the distances between central and lat- 
eral spheres. This is inspired by a 3-bead linear 
chain theoretical model [5j which has never been 
realized experimentally 



2 Materials and Methods 

2.1 Optical tweezers 

The optical tweezers setup used in this work 
consists of a laser (IPG Photonics, PYL-1-1064- 
LP, A=1064nm, P maa; =l.lW) focused through 
a water immersion objective (Zeiss, Achroplan 
IR 63x/0.90 W), trapping from below. The 
laser beam is steered via a pair of acousto-optic 
deflectors (AA Opto-Electronic, AA.DTS.XY- 
250@1064nm) controlled by custom built elec- 
tronics, allowing multiple trap generation with 
sub-nanometer position resolution. Instrument 
control and data acquisition are performed by 
custom software. The trapping potential is lo- 
cally described by a harmonic spring, and the 
trap stiffness was calibrated by measuring the 
thermal displacements of the trapped beads. 
The trapping stiffness of the central bead in this 
work is &trap = 4.1 ± 1.3 pN//im, whereas the 
external beads were held more strongly, with 

^trap, external = 8A ± 1A pN/^m. The er- 
rors quoted here are the standard deviation over 
several independent experiments, with the same 
beads and same tweezers configuration. They 
could also be a source for the slight discrepancy 
in relaxation times seen in [2](a) and (b). The 
sample is illuminated with a halogen lamp and 
is observed in bright field with a fast CMOS 
camera (Allied Vision Technologies, Marlin F- 
131B). Silica beads of 3.0/im diameter (Bangs 
Labs) were diluted to extremely low concen- 
tration to avoid any spurious bead falling into 
the laser trap. The set of three trapped beads 
was floated well above the glass slide surface (at 



about 10 times the bead diameter) to minimize 
any hydrodynamic drag from the solid surface. 
Except where specified differently, we used a so- 
lution of glycerol (Fisher, Analysis Grade) 51% 
by weight in water (Ultrapure grade, ELGA) 
which has a viscosity of 6.23 mPa s [5J. Exper- 
iments were performed at 25°C. 

2.2 Experimental protocol 

The experimental protocol is composed of two 
parts. In the first calibration stage all the traps 
are kept at rest, and the beads undergo only 
Brownian motion confined by the traps. The 
driven dynamics occurs only during the second 
stage. Data is acquired at 48.79 frames per sec- 
ond, with exposure time 19 x 10~ 3 s. A run 
lasts ten minutes, during which we collect about 
30000 frames, equally divided between calibra- 
tion and dynamics. We collected 4 runs for each 
set of parameters. Images are analysed using a 
correlation filter with a kernel optimized to the 
bead profile, followed by a 2-d least-square fit. 
This gives the center-of-mass coordinates of the 
beads in each frame with an error of the order 
of 10 _3 /^m. 

2.3 Numerical simulations 

Numerical simulations were performed by inte- 
grating directly the equation of motion of the 
three beads with Taylor's method [7J. We con- 
sidered both simulations that included thermal 
motion (i.e. adding a white noise term to the 
equations) and deterministic ones, finding no 
relevant differences. The results in the body 
of the paper mainly refer to the deterministic 
model. 



3 Experimental 

We study a swimmer composed of three beads 
immersed in a high viscosity Newtonian liq- 
uid. The beads are controlled by an opti- 
cal tweezer (8J [9] and interact with each other 
through the fluid. The laser traps are set-up fol- 
lowing the scheme of figure [T](a,b). The central 
trap is kept at rest, holding the central bead in a 



soft harmonic potential of stiffness /^rap' wrme 
the lateral ones switch between two positions. 
By forcing the lateral beads to move, we imple- 
ment the four time-reversal symmetry breaking 
phases shown in figure [Th. 

Two of the parameters that characterize the sys- 
tem are purely geometric: e, the maximum os- 
cillation amplitude for the lateral beads, and d, 
the distance between the laser positions at the 
starting phase of the cycle. Their physical in- 
terpretation is that e sets the strength of the 
stroke, while d sets the magnitude of the hydro- 
dynamic interaction. The additional parameter 
is r, the switching time of the laser, which cor- 
responds to the temporal length of each phase 
in the cycle. Finally, the viscosity rj, which is set 
in the preparation of the sample, characterizes 
the fluid, the hydrodynamic interactions, and 
the amplitude of thermal motion. The former 
parameters, together with the radius R of the 
beads, characterize completely the system. 

Despite its simplicity, the model encapsulates 
all the essential characteristics of a swimmer. 
We study the flow generation on the fluid, us- 
ing the central bead as a probe. This corre- 
sponds to studying the propulsion of a freely 
moving swimmer. By tracking the position of 
the central bead it is possible to detect a left- 
right symmetry breaking and quantify the net 
flow at varying parameters. 



3.1 Quantifying displacements 

Trapping lasers act with good approximation 
as harmonic potentials on each bead. Thus, 
the average displacement (Ax) t of the central 
bead from the position of the central trap can 
be converted using Hooke's law into the mean 
force exerted on the fluid by the whole system. 
However, a direct measurement of (Ax) t is at 
the limit of experimental resolution. The lim- 
iting factor is not the imaging resolution, but 
the presence of thermal fluctuations. This ther- 
mal noise can be compensated by longer sam- 
pling, however to reduce this error down to 
the nanometer scale would require experiments 
10 2 times longer than ours, which is not prac- 



tical Q Instead, we quantify the flow using 
a different approach. The electronically con- 
trolled trap movement temporization allows to 
precisely keep track of the phase within each cy- 
cle. We can therefore average over the repeated 
realizations, and reconstruct the mean dynamic 
cycle of the beads. It is significant to focus on 
the central one, in the stationary trap: An ex- 
ample of mean cycle for this bead is shown in fig- 
ure 121(a), where four peaks due to the lateral mo- 
tion are clearly visible. These peaks correspond 
to the maximum displacement of the central 
bead in each phase of the motion. The configu- 
rations of the three traps are indicated schemat- 
ically in the bottom panel of figure 12(a). Due 
to the different dispositions of the outer beads, 
the drags the central bead is subject during each 
phase are different and cause unequal displace- 
ments. Our method characterizes the pumping 
by means of the asymmetry of the peaks. On 
the mean cycle, we label pi, p 2 , p^, Pa (in order 
of decreasing value) the peaks in the displace- 
ment, and we define the following observables, 
devised to quantify the asymmetry of the mo- 
tion: (1) S inv := (pi - p 2 ) - G°4 - Pa), the dif- 
ference between the lower and the upper peak 
values, (2) 5 2 '■= Pi — £>4, the difference between 
the upper and the lower maximum peak val- 
ues. For an illustration we refer to the mid 
panel in figure [21(a). The reason for using both 
definitions (1) and (2) to quantify asymmetry 
is the fact that the mean cycles have different 
shapes for different values of r. For example, 
in the case of small values of r, the intermedi- 
ate peaks are not well distinguished, as shown 
in the top panel of figure [2(a). Thus, we adopt 
the observable 5 2 which makes use of the up- 
per and the lower maximum peak values only. 
In the opposite case of high values of r, see 
bottom panel of figure 121(a), the four peaks are 
well defined, and it is more effective to use the 
observable 5 inv . The advantage is that 5 inv is 
determined independently from the equilibrium 
position xq, which is needed for S 2 and cannot 
be measured during the active motion. Thus, 
it allows for more direct measurements. These 



* Already at times of the order of a few tens of minutes 
there can be sources of mechanical noise leading to drift. 



observables lack a simple physical meaning and 
should rather be thought as order parameters 
useful in characterizing the asymmetry of the 
peaks. Both these quantities change sign by re- 
verting the sequence of steps and allow to verify 
if the system behaves as expected, pumping in 
the opposite direction. We also define a third 
observable (3) s :— \p\\ + \p%\ + \ps\ + \pi\, the 
sum of absolute values of the peaks, s is al- 
ways positive, and quantifies the amplitude of 
the swimmer's motion. 

3.2 Comparison of experiments 
with simulations 

We compare the measured flow to simple nu- 
merical simulations (see methods). The hydro- 
dynamic interaction between the spheres is de- 
scribed by the Oseen tensor [TO], corresponding 
to the limit of point force, or far field. Along the 
x axis of the swimmer, the equation of motion 
for bead n is 



X , 



with 




1,2,3, 



where 7 = QirrjR is Stokes' drag, F { is the exter- 
nal force exerted on the i-th particle, and r nm in- 
dicates the relative distance between the n and 
m. The experimental values of rj and ^trap were 
used in the simulations. 

As can be seen by inspecting figure [2(b), the ob- 
served relaxation times are shorter in the exper- 
imental data of bead motions than in the sim- 
ulations at corresponding parameters. To ex- 
clude possible errors, we performed the follow- 
ing tests. First, we simulated a single particle 
in a stationary trap, including thermal fluctua- 
tions, calculated the relaxation time in the time 
autocorrelation function. For an isolated par- 
ticle, the relaxation time is the ratio ^trap/7- 
Using the experimental values in the simula- 
tions, we found excellent agreement for beads 
both in pure water and water/glycerol. Sub- 
sequently, we compared the simulation (with 
and without thermal noise) to the experiment 
for one trapped particle in water/glycerol, in 
a trap undergoing repeated displacements, see 



figure 2(c). Also in this case, the simulation 
and experiment match very well. In conclusion, 
there is excellent agreement in the single parti- 
cle regime, confirming experimental calibration 
and numerical methods. Therefore the slight 
discrepancy in the three-bead data must be due 
to the excessively simple theoretical description 
of multiple interacting bodies [j. In other words, 
the flow induced by each bead is not perfectly 
represented by a linear superposition of single- 
particle Oseen propagators [10]. This is not sur- 
prising, considering that this is a long-distance, 
or far-field approximation. A common proce- 
dure to correct for this effect makes use of a 
perturbative expansion in the parameter R/d. 
We checked this by implementing the first and 
second perturbative corrections [T3J [H] in our 
simulations, finding very similar quantitative re- 
sults as with the Oseen tensor. The reason for 
this is that the parameter R/d is not small, be- 
ing of the order of 1/3, so that the real dy- 
namics is not accessible perturbatively q De- 
spite this limit, the simulations do allow us to 
compare the experimental results with a predic- 
tion for the net force, or flow, generated in the 
fluid by the swimmer. A mapping between each 
of the observables (<5 2 , 5 inv ) and the temporal 
average {Ax)t is obtained by means of simula- 
tions. Here, sampling problems are not present 
and (Ax)t is quantifiable directly by time av- 
eraging with arbitrary accuracy; 6 2 and 5 inv are 
obtained from the steady-state mean cycle, sim- 
ilarly to experiments. For both 5 2 and 5i nv , we 
verify that there exists a one to one mapping 
with {Ax) t - We use this relation to associate as 
a function of e each point of the experimental 
curve to the corresponding mean displacement, 
obtaining {Ax) t as a function of e that is con- 
verted into a mean force via Hooke's law. The 



^Note that we can exclude that inertia and the prop- 
agation time of hydrodynamic interaction plays a role, 
as deviations were found to be at the nano-second 
scale [HI [12], well outside from the time scales of our 
experiments. 

* While approaches that do not make use of the 
explicit form of the solution of hydrodynamic equa- 
tions [151 116] may help solving this problem, in general, 
to compensate for this small discrepancy one would have 
to solve Stokes' equation with the proper boundary con- 
ditions p"7] . 



procedure is illustrated in Figure SI of supple- 
mentary methods. 



4 Discussion 

Our first result is that the system is able to 
generate flow. This is qualitatively visible in 
the distribution of position of the central bead 
(figure [Qc)). This effective potential is com- 
posed of two contributions: the harmonic po- 
tential exerted by the tweezer, and the configu- 
rational bias induced by the interaction with the 
lateral moving beads. The latter contribution 
causes the emergence of two effective asymmet- 
ric minima. This asymmetry can be quantified 
accurately using the above-defined observables 
(figure [3]). Simulations show that there exists 
a one-to-one map between S inv and <5 2 , and the 
values for the mean force, so that it is possible 
to convert 5i nv and 62 directly to forces. For 
fixed values of r, we ran a series of experiments 
to characterize the dependence of the mean po- 
sition from the displacements e of the lateral 
beads. Figure[]Jb) shows one of the two possible 
time-reversal breaking sequences. By reverting 
the sequence into the specular one, the swim- 
mer should move in the opposite direction. Ex- 
perimentally, we implemented both sequences, 
and found that this property holds. Figure [3|a) 
shows a comparison between the experimental 
and simulated observable Si nv as a function of 
e. The absolute value of the corresponding 
mean force is shown in figure E](b). The max- 
imum mean forces reached in our experiments 
are of the order of 0.03pN, roughly correspond- 
ing to a net swimming/pumping speed of about 
0.2/im/s. Note that this force is exceedingly 
small to be measured with an optical tweezer 
using conventional techniques, which has posed 
a barrier in past investigations [18]. The indi- 
rect technique used here enables to overcome 
this barrier, and can possibly be useful in dif- 
ferent contexts. 

It is interesting to study the influence of r, 
by fixing the value of e, and in figure [3]^c) we 
plot the propulsion force at varying r, compared 
with simulations. In both cases there is a good 
quantitative agreement between the model with 



Oseen interactions and the experimental data. 
By increasing the distance d between the beads, 
the asymmetry in the peaks becomes more dif- 
ficult to detect. However it is still possible 
to characterize the dependence of fluid pump- 
ing on the parameter d by using the amplitude 
variable s. In figure [3](d), s is plotted versus 
the distance, for both simulations and data. In 
both cases this quantity decays as 1/d. Quan- 
titatively, there is a small systematic deviation 
between experiment and simulations. The 1/d 
decay directly depends on hydrodynamic inter- 
actions, and can be understood with a simple 
argument on the maximum deviation of an ini- 
tially resting trapped bead subject to the per- 
turbation induced through the fluid by another 
bead at distance d relaxing for a stretch e in 
another harmonic potential (see Supplementary 
Material for further details). 



4.1 Comparison with 
model systems 



other 



It is instructive to compare our swimmer with 
the closely related theoretical model originally 
proposed in the literature [5] by Najafi and 
Golestanian, and the variants that have been re- 
cently explored [TniED]. There are two main fea- 
tures characterizing our model. Firstly, beads 
are subject to the compliance of the trapping 
potentials. The original Najafi-Golestanian 
swimmer has rigid links. More recently, the 
model has been solved for the case of chemi- 
cal transitions between states [19] and for gen- 
eral bead and link sizes and imposed deforma- 
tions [T7J [20]. However, an actuated move- 
ment by quadratic potentials has not been ad- 
dressed explicitly in the previous literature. 
While the behavior of the simplest rigid Najafi- 
Golestanian swimmer is essentially character- 
ized by the ratio e/d [5j [21], more general 
models can have more complex dynamics, im- 
posed by additional relevant time or length- 
scales. Our system is also subject to an ad- 
ditional parameter, the ratio — , with r : = 
7/^trap °f the imposed displacement time to 
the characteristic relaxation time of the bead. 
Hence the phenomenology of our system de- 
pends also on the ratio of the two time scales. 



The most important consequence is that our 
simulations show a power-law dependence of 
the speed, or mean force, on e, with an expo- 
nent that increases monotonically with increas- 
ing r, slightly larger than the quadratic law 
found for the Golestanian swimmer. Experi- 
mentally, it is problematic to measure incon- 
trovertibly this scaling, due to the large errors 
in the fit for the exponents, whose average val- 
ues, however, do exceed two. In the limiting 
case of small r, the two models behave in the 
same way. This can be understood in the fol- 
lowing way. Assume for simplicity that the lat- 
eral forcing beads can be described as free par- 
ticles relaxing trough the minimum of an har- 
monic potential. Their position follow a simple 
temporal law governed by an exponential de- 
cay x(t) ~ xo exp(— t/r ), and accordingly the 
velocity v(t) ~ (xq/to) exp(— t/r ). For small 
values of r, the ratio t/tq is small ensuring that 
the velocity of the trapped bead does not de- 
cay significantly during each step in the cycle, 
being a constant proportional to tq as in the 
Najafi Golestanian swimmer, where the beads 
have a constant displacement velocity [5j. On 
the other hand, our results confirm that in the 
small t limit the mean force scaling with e/d is 
compatible with a quadratic law. This regime 
is also the one where the swimmer's propulsion 
is optimal. In the numerical experiments, the 
fitted exponents we find are consistently larger 
than 2 (between 2.12 and 2.41). However, the 
predicted trend with e cannot be confirmed ex- 
perimentally because of the large errors (Sup- 
plementary Figure S2). The second difference 
with the original Najafi-Golestanian swimmer 
is that our swimmer is not intrinsic, but actu- 
ated by external forces [U |22]. In particular, 
the total actuating force is not null instanta- 
neously, but only over a cycle. This has been 
shown to lead to qualitatively distinct behavior 
in other systems [23]. In our case, using sim- 
ulations of an analogous intrinsic swimmer ac- 
tuated by two-body two-state springs, we find 
that the differences are mostly quantitative, the 
intrinsic variant being slightly faster and more 
efficient (Supplementary Figure S3). 

Finally, it may be of interest to compare this 
pump to other engines. We can estimate the 



dissipation rate in our pump to be between 
0.5 x 10~ 17 and 4 x lCT 17 J/s, which is com- 
parable to an E. coli bacterium swimming in 
water [23]. The efficiency of this pump is of 
the order of 10~ 3 % (calculated as the ratio 
of the mean output to input power) which is 
lower than for typical biological systems such 
as E. coli pTJ, [25]. These results are important 
on fundamental grounds and for applications in 
nano-scale machines. 
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Fig. 1 Experimental realization of the low Reynolds number swimmer, and proof 
of force generation, (a) Scheme of the lasers' disposition during a single phase of 
motion, in which the active traps (green) act on the beads as harmonic potentials. 
The parameters e and d are sketched, (b) Sequence of snapshots of the experiment 
showing four time-symmetry breaking steps in the basic cycle. Crosses indicate the 
active laser trap positions. By moving the lateral beads and keeping the central one 
at rest, it is possible to study the flow generated on the fluid, using the central bead 
as a probe. A left-right asymmetry in its position signals that flow is generated. Top 
panel: starting configuration, where the distances between beads take their maximum 
values, (c) Two-dimensional map of the effective potential felt by the central bead. 
The plot is obtained by taking the logarithm of the relative frequencies of the positions 
occupied by the central bead during 500 repetitions of the basic cycle. The colour scale 
indicates in dark the most frequent position. The map shows a left-right symmetry 
breaking (the left-hand minimum is deeper) , which is a proof of the mean generated 
flow in the system. 

Fig. 2 Reconstruction of the mean cycle. The experiment is a sequence of repetitions of 
the elementary cycle illustrated in figure [TJ By treating each period as a single realiza- 
tion and averaging over all realizations, one recovers the mean cycle of the dynamics, 
averaged over the thermal fluctuations due to Brownian motion in the fluid, (a) Detail 
of the mean cycle for the central bead (<^ with errorbars) compared with simulations 
(solid lines) for different r. The plot reports the longitudinal displacement coordinate 
of each bead as a function of the cycle phase in the interval (0, 2tt). The peaks become 
more defined as r increases. Top panel: r = 40 ms ; middle panel: r = 80 ms ; bottom 
panel: r = 320 ms. The position of the four peaks and the traps' configuration in the 
four phases of the cycle are illustrated in the middle and bottom panel respectively. 

(b) Mean cycle for the whole system (three beads). The longitudinal displacement 
coordinate of each bead is plotted as a function of time. Dotted lines correspond to 
experimental values for the left bead (black), the central bead (red) and the right bead 
(green). These are compared with the results of numerical simulations (solid lines with 
the same color code), (c) Close up and complete cycle showing relaxation of a single 
bead in a switching trap, with no interacting partners. 

Fig. 3 Quantitative characterisation of the swimmer. The main observable is the 
mean force felt by the central bead, measured by its mean displacement in the optical 
trap. The mean force is studied here as a function of the parameters e and r. (a) 
shows Si nv as a function of the lateral amplitude e, comparing experimental values 
(triangles) with simulation (continuous line). There are two possible time-reversal- 
symmetry breaking cycles, corresponding to the two opposite pumping directions. By 
temporally reversing the cycle showed in figure (1), 5i nv changes sign, showing that 
the direction of pumping can be controlled. The following panels show experimental 
data (<0) compared to simulations (continuous line). Each <0> represents the average 
over 4 different experiments and the error bars correspond to the standard deviations 
of the four values, (b) Shows the mean force (obtained from the observable 5 inv and 
numerical simulations) as a function of e, for d = 6/im, r = 80 ms, rj = 6.23 mPas. 

(c) Mean force as a function of r. Here the other experimental parameters are d = 6/im, 
e = l//m, 7] = 6.23 mPas. (d) Amplitude of motion as a function of d. Increasing d 
the force becomes weaker and difficult to detect, and we have to use the amplitude 
parameter s to compare with simulations. The measured and simulated values of s 
scale as d" 1 (as discussed further in the Supplementary Material). 
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Figure 1: Experimental realization of the low Reynolds number swimmer, and proof 
of force generation. 
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Figure 2: Reconstruction of the mean cycle. 
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Figure 3: Quantitative characterisation of the swimmer. 



SUPPLEMENTARY MATERIAL 

Additional considerations on the model and the data anal- 
ysis 



Scaling with d 

We present here a simple argument for the scaling of the amplitude s with distance d, explaining 
the scaling seen for experiments and simulations in figure 3(d). 

To estimate the (^-dependency of the amplitude during a cycle, it is sufficient to describe the 
displacement of the central bead in one of the four sub-phases. We consider the simplified but 
physically equivalent situation of two interacting trapped beads at distance d. Initially the left 
bead, which represents the central bead in the experiment, is at rest and the right bead is out of 
its equilibrium position by a distance e, as it would be if in the instant the trap has been shifted. 
The initial conditions are thus xl(0) = 0, Xr(0) = d and the equilibrium positions of the trapping 
potentials are Xo,l = 0; rr 0i _R = (d — e). To a first approximation, ("zeroth order" in R/d) the left 
bead is still, and the position of the right bead follows a simple relaxation law 

4 0) W= e (e---l)+rf • (2) 

This unperturbed solution can be used to estimate (by the force balance with the fluid) the source 
of force F R = — fctrap( a; it — °0 applied by the right bead on the fluid during its relaxation. This 
can then be used into equation (1) of the main text for the left bead, giving 

To i-L = x l ~ ~o1 Xr ' ft' 

where x R appears as an external perturbation, and we approximated the distance between the 
beads with d, which is justified in the limit of large distances d » e. With this assumption, the 
problem becomes linear and the solution to order R/d can be calculated as 

m , , 3 -Re t __L , ,, 

*■■ (f) = mV ' ■ (4) 

In turn, this solution could be used as a source for the equation for xr, to obtain hierarchically 
the higher order contribution in R/d to its motion. The value of the peak of the central bead in 
the experiment s can be estimated by the maximum displacement of the left bead 

T (i) _ ^3l (K) 

This argument implies that the leading order scaling of each peak in a cycle, and hence of s, is 
1/d. The argument has the advantage of showing how the hydrodynamic interaction tensor comes 
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into play explicitly after a trap switches its position. 

In the same linear approximation, the coupled equations (eq.l) from the main text can even be 

solved directly in a straightforward way, and the resulting maximum displacement is: 

XL = ^! I (6) 

which has the same behavior in the limit of large d/R. 

Strictly speaking, this result is applicable in the regime r >> r , in which the beads have the 
time to fully relax in the trap potentials. In the opposite limit, r << r , since the subcycle ends 
while the central bead is still moving away from the center of its trap, the maximum position can 
be estimated using the same solution, x L (t), by the position assumed by the central bead at the 
end of this subcycle, i.e. the instant t — r. Thus as 

/(1) _ 3&T i 
Xl ^~ 2dr 6 ' {7) 

and the scaling with d is unaffected. 
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Relating the observables 62 and Sinv with the mean force. 

We summarise here how the procedure to relate the experimental data to the induced flow. 
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conversion into mean force 
(Hooke's law) 




Supplementary Figure SI: Converting observables into mean force. This scheme illustrates 
the procedure used to obtain the mean force from the experiments, by means of simulations. 
Step 1: Because of the different shapes of the mean cycles for different values of r, we define two 
different observables 62 and Sinv to quantify the asymmetry in the displacements for the central 
bead. The peaks which enter in the definitions are indicated with red arrows. Step 2: (a) We 
analyze the mean cycle of experiments and simulations, extracting each observable at varying e. 
We also calculate directly the temporal average (Ax) for the position of the central bead, (b) 
Comparing these results and eliminating the dependence from e, we find that there exist a one to 
one mapping between each observable and the temporal average (Ax) . (c) Using this relation the 
curve of each observable as a function of e can be converted into a curve of the mean position as 
a function of e. Step 3: using Hooke's law we convert the mean displacement into a mean force 
as a function of e. 
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Scaling with r 
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Supplementary Figure S2: Scaling law for the mean force at varying r. (a) Plot on log- log 
scale of the simulated mean displacement from the equilibrium position of the central bead. The 
curve is a power law with different exponents for different values of r. In our model the exponent 
varies, and increases monotonically for increasing r. In the limit of small r the mean displacement 
follows a power law with exponent close to 2 (dashed line), which resembles the behavior of the 
Golestanian swimmer, (b) Comparison between experimental and theoretical data for r = 320 
ms on log-log scale. Due to the large error bars on the experimental curve, the determination of 
the exponent from the experiment is subject to large uncertainties. 
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Intrinsic Swimmer 



Figure [S3] shows a comparison of our model with simulations of an analogous intrinsic swimmer. 
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Supplementary Figure S3: Comparison of the propulsive forces for extrinsic and intrinsic 
swimmers. The extrinsic swimmer is studied experimentally and numerically in this work, 
whereas we can only study the intrinsic swimmer, actuated by two-state springs, numerically. 
The plot shows the mean force with varying e, for simulations of the two models, with parameters 
d = 6fim t = 80ms. The difference between the two swimmers is only quantitative. 
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